Transverse momentum dependence of the quark helicity distributions and the Cahn 
effect in double-spin asymmetry All i n Semi Inclusive DIS 



> 
O 



M. Anselmino 1 , A. Efremov 2 , A. Kotzinian 2,3 - 4 and B. Parsamyan 3 

1 Dipartimento di Fisica Teorica, Universitd di Torino, and 
INFN, Sezione di Torino, Via P. Giuria 1, F10125 Torino, Italy 

2 JINR, 141980 Dubna, Russia 

3 Dipartimento di Fisica Generale, Universitd di Torino, and 
INFN, Sezione di Torino, Via P. Giuria 1, 1-10125 Torino, Italy 

4 Yerevan Physics Institute, 375036 Yerevan, Armenia 

\£l • Within the LO QCD parton model of Semi Inclusive Deep Inelastic Scattering, IN — » ihX, 

with unintegrated quark distribution and fragmentation functions, we study the Put dependence of 
f*"^ , the double longitudinal-spin asymmetry All- We include 1/Q kinematic corrections, which induce 

£NJ . an azimuthal modulation of the asymmetry, analogous to the Cahn effect in unpolarized SIDIS. 

We show that a study of All and of the weighted DSA A c L ° L "*' h allows to extract the transverse 
momentum dependence of the unintegrated helicity distribution function g\ L (x, k±) [or Aq(x, k±)]. 
y-^ . Predictions, based on some models for the unknown functions, are given for ongoing COMPASS, 

■ HERMES and JLab experiments. 
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PACS numbers: 13.88.+e, 13.60.-r, 13.87.Fh, 13.85.Ni 
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00 . I. INTRODUCTION 

-xf ; 

qq ■ After many years of intensive experimental and theoretical study the partonic origin of the nucleon spin still remains 
| mysterious. In particular, while a lot of understanding has been achieved concerning the longitudinal structure of a 
\& • fast moving proton - the ^-dependence of the parton distribution functions and of the helicity distributions - very 
little is known about the transverse structure. Transverse refers to the direction of motion and concerns both the 
"P^ , transverse spin distributions and the parton intrinsic motion, k±. These unexplored degrees of freedom cannot be 
O 1 considered as minor details in our modeling of the nucleon. Without a good knowledge of the total intrinsic momentum 
carried by the partons, and its connection with the spin, we could never understand the parton orbital motion and 
J^" 1, progress towards a more structured picture which goes beyond the simple collinear partonic representation. 
(— i | Recent data on single spin azimuthal asymmetries in semi-inclusive deep inelastic scattering (SIDIS, IN — > ihX) 
obtained by HERMES Q, COMPASS and CLAS-JLab 0] collaborations triggered a lot of interest towards the 
. £h ' transverse momentum dependent (TMD) and spin dependent distribution and fragmentation functions (PDFs and 
FFs). In fact spin-fe^ correlations induce new spin effects, which would be zero in the absence of intrinsic motion and 
' are related to transverse polarizations. On the other hand, it has been known for a long time that the dependence on 
the transverse final hadron momentum, PhT, observed in unpolarized SIDIS processes, can be related, for PhT = \Pht\ 
values up to about 1 GcV/c, to quark intrinsic motion 0, @, Q. 

We consider here polarized SIDIS processes, at twist-two in the parton model, with transverse momentum dependent 
distribution and fragmentation functions. Such processes can be described in terms of six time reversal even 0, HJ 
and two (naively) time reversal odd PDFs. The dependence on partonic intrinsic motion induces a dependence on 
PhT- bi addition, at 0(k±/Q), kinematic corrections induce a dependence of the unpolarized cross section on the 
azimuthal angle <ph between the leptonic and the hadron production planes - the so called Cahn effect Q. It was 
shown in Ref. Q that a careful study of the dependence of the cross section on the final hadron momentum allows to 
extract the average values of intrinsic momenta in unpolarized PDFs and FFs. 

We expand on the work of Ref. Q and evaluate the role of partonic intrinsic motion in polarized SIDIS; in particular, 
on the double spin asymmetry (DSA) for the scattering of longitudinally polarized leptons off a longitudinally polarized 
target, All, where longitudinal refers to the incoming lepton direction, in the laboratory frame. We show that a 
study of All and of the weighted asymmetry A L °^ h allows to learn about the k± dependence of the quark helicit 
distribution g^ L (x, fcj_) [or Aq(x, k±)]. Similar results for the TMD DF g^ T have been obtained in a recent paper 

by considering the longitudinal-transverse DSA A L ° T ^ h ~^ s \ 

The article is organized as follows. In Sec [IT] we shortly recall the relevant formalism for polarized SIDIS. In Scc lIIII 
some predictions for the double longitudinal spin asymmetries are presented. The results are given for different sets of 
kinematical cuts, according to the experimental setups of HERMES, COMPASS and JLab experiments; they indicate 
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the best kinematical regions for the asymmetry to be sizeable. Finally, in Sec. II VI we shortly discuss our results and 
draw some conclusions. 

II. POLARIZED CROSS SECTION 

Following Ref. 0, we consider the polarized SIDIS in the simple quark-parton model, with unintcgrated parton 
distributions. We use the standard notations for DIS variables: I and I' are, respectively, the four-momenta of the 
initial and the final state leptons; q = I — £' is the exchanged virtual photon momentum; P (M) is the target nucleon 
momentum (mass), S its polarization; P] % is the final hadron momentum; Q 2 = —q 2 ; x = Q 2 /2P ■ q; y = P ■ q/P ■ £; 
z = P ■ Pn/P ■ q, Q 2 — xy(s — M 2 ). s = (£ + P) 2 . We work in a frame with the z-axis along the virtual photon 
momentum direction and the x-axis in the lepton scattering plane, with positive direction chosen along the lepton 
transverse momentum. The produced hadron has transverse momentum Put] its azimuthal angle, <ph, and the 
azimuthal angle of the transverse nucleon spin, (f>s, are measured around the z-axis (for further details see Ref. Q). 

We consider longitudinally polarized protons and leptons, where longitudinal (according to the laboratory setup) 
refers to the initial lepton direction. It then results that a proton with longitudinal spin S along the incoming lepton 
direction, has a transverse - with respect to the 7* direction - spin component: 



5 T = 5sin6» 7 , (1) 



where 



AM 2 x 2 ( M 2 x 2 y 2 \ 2MxJY=y 
sm 6, = \\ Q2+ 4M2x2 [l-y- Q2 ) * q ( 2 ) 



This component gives contributions of order M/Q. 

Keeping only twist-two contributions and terms up to 0(M /Q) the cross section for SIDIS of longitudinally polarized 
leptons off a longitudinally polarized target can be written as: 

d 5 a^ 2a 2 

{H fl + X(S L H giL + STU giT )} , (3) 



dx dy dz d 2 PhT xy 2 s 

where the arrows indicate the direction of the lepton (— >) and target nucleon (<=) polarizations, with respect to 
the lepton momentum; A, Sl and St are the magnitudes of, respectively: the longitudinal beam polarization, the 
longitudinal and the transverse target polarization. Notice that <= stands for a nucleon with a polarization vector, 
in the laboratory frame where the nucleon is at rest, opposite to the initial lepton momentum. For a polarization 
one reverses the signs of the Sl and St terms. 

The three terms have a simple partonic interpretation: 

n h =J2 e l J ^k±f!{x,k±)iry 2 ^^-D^z,p±), (4) 



"2 ^ 2 

=5>« / d 2 k x g\ L {xM)W S —^D h 9 {z, PA _), (5) 



n aiT = -22 e l J d 2 k ± — cos ipg^ix^^) ny 2 ^ D*(z,px) , (6) 

and deserve some comments. 

• The partonic factorized structure of the above equations is supposed to hold in the large Q 2 kinematic region 
where P^t — Aqcd - fci C [13]. It neglects terms of 0(k±/Q) 2 , in which case 

p ± = P hT - zk± , 

where p^ is the intrinsic transverse momentum of the hadron h with respect to the fragmenting quark direction. 
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• The first two contributions, Eqs. ([4| and ([5]), give, respectively, the unpolarized cross section and the helicity 
asymmetry 



(Fa 



dx dy dz d 2 Pi 



hT 



2a 2 
x y 2 s 



n 



d 5 a 



d 5 a^ 



Aa 2 



h 



dx dy dz d 2 PhT dx dy dz d 2 P^t 



,2„ I 



(7) 



x y*s 

where +, — stand for helicity states. The quark intrinsic motion induces a kinematical azimuthal dependence, 
via the elementary polarized cross section 



oc — 



A A„ 



u 2 ) 



dQ 2 dip 



t 2 



(8) 



where X q denotes the quark helicity. Keeping the terms up to order of kj_/Q the Mandelstam variables for the 
non-coplanar iq — ► iq scattering are expressed as 



fcj 



1-2^1^ 



cos ip 



-Q 



-xys, 



(9) 



-xs (1 — y) 



2k, 



COS tp 



where tp is the azimuthal angle of k±, d 2 k± = k± dk± dip. Eq. (0| then gives the unpolarized Cahn effect 0], 
while Eq. (J3J) gives the corresponding effect for the polarized (helicity) cross section, both at 0(k±/Q). 

Eq. ^ contains another cos<^ dependence, of different origin. While the distribution functions f^(x,ki_) 
and gf L (x,k±) which appear in Eqs. Q and (JSJ), are just the k± dependent unpolarized and longitudinally 
polarized (helicity) PDFs, which, upon integration over d 2 k±, give the usual fi(x) [or q(x)] and g\(x) [or 
Aq(x)] distributions, the quantity 



^ cos (pgl^,(x,k±) = A/^/St 



(10) 



is related to the number of partons longitudinally polarized inside a transversely polarized proton 0, d, EH : 
it can only depend on the scalar product between the two corresponding polarization vectors, which gives the 
cos(^>s T — tp) — — costp factor explicitly shown (sec, for example, Eq. (C19) of Ref. [HI]). This distribution 
is a leading- twist one, not suppressed by (k±/Q) small factors. However, Eq. ^ will be multiplied by St, 
which is of 0(M/Q), Eqs. ([I])-©; for this reason, in Eq. ([6]) we shall not take into account the extra (k±/Q) 
kinematical terms contained in (s 2 — u 2 ) of Eq. (J8j) . 

The integrals in Eqs. Q-© can be analytically performed, if one assumes a simple factorized and gaussian 
behaviour of the involved TMD PDFs and FFs 



gl T (x,k±) 
gl L (x,k±) 



- fi 



1 



fi( x ) — 2 ex P 
ir[i 2 



D h q {z) 
9 q i T ( x ) 



— 5- exp 

7T/4 



— * exp 

TV fl( 



— 5- exp 
TM2 



k\ 
2 
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yielding, at 0(P hT /Q): 



1 + (1 - yf - 4(2 - y) ^Y^y 



Z $ P hT 



COS (f) h 



exp 



( /In 



Hd + z 2 14 



(11) 

(12) 
(13) 
(14) 

J2e 2 q f!(x)D h q (z), (15) 



K giL = y 



2-2/-4VT - 



' (4 P hT 



Q(j** d + z*i4) 



COS 4>h 



V 2 D 



+ z 



2 nl 



E 



e 2 gt{x)D h q {z), 



(16) 
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"ftgiT = -2/( 2 - v) 



Z Hi PhT 



cxp 



M(ft 



V?) 



27 COS (j) h 



Y / e 2 q9 q 1T (x)D h q (z) 



(17) 



III. PREDICTIONS FOR All 



We use Eqs. and ([15]) — (|17p to compute observables which depend on partonic intrinsic motions. Notice that we 
have allowed different average values of (k 2 ^) for the different distribution functions: (k 2 ) = /ig for the unpolarized 
distributions, (k^_) = /i 2 , for the helicity distributions, and (k 2 ± ) = fj, 2 for gf^(x,k±); each of these value is taken to 
be constant and flavour independent. For the fragmentation functions we have (p 2 ^ = n 2 D - Following Ref. [|| we use 

Hi = 0.25 (GeV/c) 2 p? D = 0.20 (GeV/c) 2 , (18) 

while wc consider fj, 2 and fi 2 as free parameters, which can give interesting information on the quark transverse motion 
in polarized protons; the naive positivity bounds imply that we should have 

lA < Mo A < Mo • (19) 

Our approach is supposed to hold up to PhT — 1 (GeV/c) 0. Above that higher order pQCD corrections must 
be taken into account, and lead to tiny variations of the values given in Eq. (|18[) @; however, we shall consider 
experiments which are expected to produce data mainly in the low PhT region, and both our approach and /ig D values 
are well adequate. 

We consider the PhT dependence of the double longitudinal spin asymmetry 

dd>h \ da^ — da^ 1 , 
A LL {x,y,z,P hT ) = ^ — (20) 



and the coscfth weighted asymmetry, defined as 



XS J 2 ^ d(f>h [ da^ + da 



2 f 2 * d<ph [ da^ — dcT^ } cos <f>h 



A c °l^{x,y,z,P hT )= - J0 = >^ZZIL . (21) 



XS J * dtfih [ da^ + da' 



From Eqs. ([T3D-([T7J) one has 



A LL (x,y,z,P hT ) = ^±, (22) 

^0 



with 



and 



Aa LL = V{2 2 V) 2 1 2 2 exp (- 2 We 2 q 9l{x)D h q {z) . (23) 



Analogously, Eqs. (JT5l)-([T71) and © give 



A?^(x,y,z,P hT ) = A(Jll +A(Jlt , (25) 



CO 

where the contribution from the longitudinal part of the target polarization is given by 

,2 v. _ / r>2 



A(T cos^ = _ 4 ^L^ z ^^t (_ P hT We 2 gi (x)^(z), (26) 

xy Q(^ 2 D + z 2 ^ 2 ) 2 l \ ^ + 
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and the contribution from the transverse part of the target polarization by 

Aa cos, h = -2{2-y)^T—y z»\P hT ( P 2 hT \y, . 

iT 2/ Q(ji 2 D +z 2 fi 2 ) 2 V + -z Mi / g 9 

Of course, both the numerator and denominator of Eqs. ((20)) and (|2Tj) can be integrated over some of the variables, 
according to the range covered by the setups of the experiments we shall consider: 



COMPASS: positive (h+), all {h) and negative (h~) hadron production, Q 2 > 1.0 (GeV/c) 2 , W 2 > 25 GeV 2 , 
0.1 < x < 0.6, 0.5 < y < 0.9 and 0.4 < z < 0.9 



• HERMES: tt+, tt° and tt" production, Q 2 > 1.0 (GeV/c) 2 , VP 2 > 10 GeV 2 , 0.1 < x < 0.6, 0.45 < y < 0.85 and 
0.4 < z < 0.7 

• JLab at 6 GeV: tt+, tt and tt~ production, Q 2 > 1.0 (GeV/c) 2 , W 2 > 4 GeV 2 , 0.2 < a: < 0.6, 0.4 < y < 0.85 
and 0.4 < z < 0.7. 

We start by considering Eqs. (|2"2")) - (f2"4")) . Notice that they are leading-twist quantities, not suppressed by any inverse 
power of Q. Concerning the usual integrated distribution and fragmentation functions we use the LO GRV98 fl2T ] 
unpolarized and the corresponding GRSV2000 [l3[ polarized (standard scenario) DFs, and Kretzer [l4| FFs. We can 
then compute the Pkt dependence of All, depending on the only unknown quantity We plot the results of our 
computations in Figs. 1 and 2, for a proton and dcutcron (+ neutron, for JLab) target, respectively. 
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FIG. 1: Predicted dependence of All on PhT, for scattering off a proton target, with different choices of 0.1 (GeV/c)^ 
continuous, 0.17 (GeV/c) 2 - dashed and 0.25 (GeV/c) 2 - dot-dashed lines. 



The results depend clearly on the relative values of (k 2 ^) for the unpolarized and helicity distribution, /Xg and 
respectively: All(Piit) is approximately constant if /x 2 , = [Iq = 0.25 (GeV/c) 2 , whereas it sharply decreases with 
PhT if A*l < A*o- T ne trend of All^Pht) is thus a significant indication of the average quark transverse motion inside 
unpolarized versus longitudinally polarized nucleons. Although our numerical estimates are based on the gaussian 
factorization ansatz, Eqs. (11)-(14), we expect them to have a more general interpretation and information content. 
The PhT dependence of All reflects, essentially, the difference between the k± dependence of ff{x, k±) and g\ L { x i 
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FIG. 2: Predicted dependence of All on Put, for scattering off a deuteron (and neutron for JLab) target, with different 
choices of 0.1 (GeV/c) 2 - continuous, 0.17 (GeV/c) 2 - dashed and 0.25 (GeV/c) 2 - dot-dashed lines. 



independently of their functional forms; the trend of All(Piit), whether constant or decreasing, reveals the behavior 
of g\ L {x, fcj_)/ fi(x,k±) and their relative k± dependence. 

Similarly, we can use Eqs. (|M)) - P7|) in order to give some estimates of A c L °^ h . Notice that Aa L °l'^ h and Act^t^ are 
(kinematical) higher-twist quantities, proportional to Pht/Q', in addition, contains one unknown function, 

namely gf T (x), related to the hclicity distribution of partons inside a transversely polarized proton. In the absence 
of any better guidance, we adopt the same strategy as in Ref. [||. We start by noticing that, from Eq. (fTF 



gl±"(x) = d<k ± -^g$(x,>4) = ift.] . (28) 



According to Refs. [1, [l5| g^\x) is directly related to the DF g%(x), which has both twist-two and higher- twist 
contributions, 

9l{x) = ^g q ^\x). (29) 

This relation, although much debated, arises from constraints imposed by Lorentz invariance on the antiquark- 
target forward scattering amplitude and the use of QCD equations of motion for quark fields Q. If, in addition, one 
uses the Wandzura and Wilczek (l6j approximation for the twist- two part of g^ix), 

gl{x)^-g\{x)+ f'dx' 9 -^!, (30) 



the following relation can be derived [17j . 



which, via Eq. (|28p . allows to express g1 T (x) through the well known integrated hclicity distributions. 



g$\x)*x fdx'i^l, (31) 
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Although such a procedure is appealing and convenient, wc should stress there are strong arguments [l8|-[20( 
against the validity of the relation (f2T))) . Therefore, we should consider the above expression, Eq. (|3~lj) , only as a rough 
model for the otherwise unknown function g\ T (x). 
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FIG. 3: Predicted dependence of A C j°^ h on P^t for scattering off a proton target with different choices of 0.1 (GeV/c) 2 - 
continuous, 0.17 (GeV/c) 2 - dashed and 0.25 (GeV/c) 2 - dot-dashed lines. Each line splits into three almost overlapping lines 
corresponding, for each value of to three different values of ^\ = (up-down) 0.1, 0.15 and 0.2 (GeV/c) 2 . 



In Fig. [3] we show our predictions for A c £? (Pht) as measurable by COMPASS, HERMES and JLab collaboration 
experiments on a proton target. The analogous results, for scattering off a deuteron target (and a neutron target as 
well, for JLab) are shown in Fig. [4j Again, we present the results for three different choices of /i 2 , = 0.1, 0.17 and 
0.25 (GeV/c) 2 , which turn out to be well different from each other. Instead, when varying the values of /xf our results 
hardly change: each line, obtained at a fixed /i 2 . value, simply splits in three almost overlapping lines (corresponding, 
from up down, to fif = 0.1, 0.15 and 0.2 (GeV/c) 2 ). This is not surprising, as, when adopting the expression (1281) . 
there remains little dependence on \j\ in Eq. (|27p . Our computations show instead a clear strong dependence on /i 2 ,. 

It is interesting also to compute the dependence of A c £^ h on each of the other single variables; for example, the 
x-dependence is computed as 

/;r— dPl T Jdyfdz(Aal°^ + A*%+») 

AT L * h (x) = hT — p2 , (32) 

dP hTldyJdza 

while the y- and z-dependences are calculated in a similar way. In Figs. [5] and [7] we present the x-, y- and z- 
dependences of A c °l <t ' h integrated over P hT with P hT , min = 0.5 GeV/c and /i 2 =0.15 (GeV/c) 2 , /i| =0.25 (GeV/c) 2 
for COMPASS {P hT , ma x = 2 GeV/c), HERMES (P hT , max = 1-5 GeV/c) and JLab (P hT , ma x = 1 GeV/c) kinematics. 



IV. DISCUSSION AND CONCLUSIONS 



We have studied the P hT dependence of A LL and A™^ 0h , measurable in SIDIS processes by COMPASS, HERMES 
and JLab collaborations. For P^t values up to ~ 1 GeV/c this dependence is entirely generated by intrinsic motion, 
both of partons inside the nucleons, and of hadrons in the parton fragmentation process 0, 
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FIG. 4: Predicted dependence of A c ^^ h on Put for scattering off a deuteron (and neutron for JLab) target with different 
choices of 0.1 (GeV/c) 2 - continuous, 0.17 (GeV/c) 2 - dashed and 0.25 (GeV/c) 2 - dot-dashed lines. Each line splits into 
three almost overlapping lines corresponding, for each value of (x 2 ., to three different values of fjif = (up-down) 0.1, 0.15 and 0.2 
(GeV/c) 2 . 



Within a simple factorized gaussian model for the kj_ and p± dependence of the distribution and fragmentation 
functions, it turns out that All^P^t) is strongly sensitive to the relative value of in unpolarized (xXq) and 
helicity (/x^) quark distributions: similar values, /Xq ~ /i^, would reflect into an approximately constant A^^Phr), 
while ^2 < /Uq' WOIU d lead to a decreasing trend. Such different behaviours are expected in general, independently of 
the factorized gaussian assumption, as the shape of A^^Phr) is essentially related to the ratio of the k± dependence 
of g\ and j\. Notice, however, that we have assumed the same constant values of (k\) and for all quark flavors; 
more involved choices might lead to different behaviours. A comparison of the quark intrinsic transverse momentum 
in unpolarized and longitudinally polarized protons might give new important information concerning the spin and 
orbital motion of quarks. For example, one expects that parton transverse motion contributes to the longitudinal 
component of the angular momentum, differently inside unpolarized and longitudinally polarized nucleons. 

The Pht dependence of A c ^^ h is not only related to kinematical non-collinear contributions, but also to a new TMD 
and spin dependent function, which gives the number density of longitudinally polarized quarks inside a transversely 
polarized nucleon. This function, gf^, induces a cos</>/j dependence, but it is unknown; we adopted a much debated 
relationship, together with the twist-two part the Wandzura-Wilczck sum rule (and the usual Gaussian factorization), 
in order to link the x-dependent part of gf^ to the integrated helicity distributions. Within such an approach, it 
turns out that also A c ^^ h (P^t) has a strong dependence on alone, thus giving further information on the average 
transverse motion of quarks inside a longitudinally polarized proton. 

We conclude by noticing, as it was done in Ref. Q, that the exact k± dependence of the distribution functions 
ff(x,k±), gf L (x,k±) and gf T '{x,k±) is crucial when considering the general positivity bounds of Ref. [2l[, which 
involve in one inequality the three previous functions and the Sivers function. The k± dependence might play an 
essential role in fulfilling the inequality, and a check of its validity is a fundamental test for the self consistency of the 
LO QCD description of SIDIS processes. 
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FIG. 5: Predicted dependence of A c ^l on x i V an d z, for proton and deuteron targets, for COMPASS. 
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FIG. 6: Predicted dependence of A^J on x, y and z, for proton and deuteron targets, for HERMES. 
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